Interparticle correlations in the production of J /ip pairs in proton-proton collisions 
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We focus on the problem of disentangling the single (SPS) and double (DPS) parton scattering 
modes in the production of J/ip pairs at the LHC conditions. Our analysis is based on comparing the 
shapes of the differential cross sections and on studying their behavior under imposing kinematical 
cuts. On the SPS side, we consider the leading-order 0{a'^) contribution with radiative corrections 
(taken into account in the framework of the fct-factorization approach) and the subleading 0{a^) 
contribution from pseudo-diffractive gluon-gluon scattering represented by one gluon exchange and 
two gluon exchange mechanisms. We come to the conclusion that disentangling the SPS and DPS 
modes is rather difficult on the basis of azimuthal correlations, while the rapidity difference looks 
more promising, provided the acceptance of the experimental detectors has enough rapidity range. 
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I. INTRODUCTION 

Since it was first observed, charmonium production in 
hadronic collisions has been a subject of considerable the- 
oretical interest. Production rates and their dependence 
on the different kinematic variables provide important 
tests for comparing theoretical models. 

In the last years, the production of J/tp pairs has at- 
tracted a significant renewal attention in the context of 
searches for double parton scattering processes A 
number of discussions has been stimulated by the recent 
measurement [2] of the double J/ip production cross sec- 
tion at the LHCb experiment at CERN. Theoretical esti- 
mates based on both coUinear [^-Q and /cj-factorization 
approaches show that the single (SPS) and double 
(DPS) parton scattering contributions are comparable in 
size and, taken together, can perfectly describe the mea- 
sured cross section. 

To disentangle the SPS and DPS mechanisms one 
needs to clearly understand the production kinematics. 
Naive expectations that the SPS mechanism should re- 
sult in the back-to-back event configuration received no 
support from the later calculations. Including the initial 
state radiation effects (either in the form of fct-dependent 



gluon distributions Q or by means of simulating the par- 
ton showers in a phenomenological way (3|) washes out 
the original azimuthal correlations, thus making the SPS 
and DPS samples very similar to each other. One cannot 
exclude, however, that the situation may change under 
imposing certain cuts on the J/^p transverse momenta. 
On the other hand, it has been suggested 0, [1] that the 
DPS production is characterized by a much larger rapid- 
ity difference between the two J/ip mesons. The domi- 
nance of the DPS contribution over SPS at large rapidity 
difference was discussed recently also for pp ccccX 
reaction The goal of the present study is to 

carefully examine the J/tp pair production properties in 
the different kinematical domains paying attention to the 
different contributing processes. On the SPS side, we 
consider the leading-order 0(af ) subprocess (with radia- 
tive corrections taken into account in the framework of 
the fct-factorization approach) and the subleading 0{a^) 
contribution from pseudo-diffractive gluon-gluon scatter- 
ing represented by one-gluon exchange and two-gluon 
exchange mechanisms; the latter mechanisms yet have 
never been discussed in the context of searches for DPS. 
On the DPS side, we consider the prompt production of 
J/ip pairs including the direct gg J/'<p9 contribution 
and radiative Xc decays. 
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II. THEORETICAL FRAMEWORK 
A. SPS contributions 

At the leading order, 0{a'^), the SPS subprocess 
g + g ^ J/ip + J /ip is represented by a set of 31 "box" 
diagrams, with some examples displayed in Fig. 1. Our 
approach is based on perturbative QCD, nonrelativistic 
bound state formalism and the /cj-factorization 

ansatz in the parton model. The advantage of us- 

ing the /ct-factorization approach comes from the ease of 
including the initial state radiation corrections that are 
efficiently taken into account in the form of the evolution 
of gluon densities. The calculation of this subprocess is 
identical to that described in Ref. Q- 

As usual, the production amplitudes contain spin and 
color projection operators that guarantee the proper 
quantum numbers of the final state mesons. Then, the 
J/^ formation probability is determined by the radial 
wave function at the origin of coordinate space |7?.(0)p; 
the latter are known from the J /ip leptonic decay width 
[l7| . Only the color singlet channels are taken into con- 
sideration in the present study since this approach was 
found to be fully sufficient [1^ to describe all of the 
known LHC data on J /ij) production. 

The evaluation of Feynman diagrams is straightfor- 
ward and follows the standard QCD rules, with one reser- 
vation: in accordance with the fct-factorization prescrip- 
tion 14], the initial gluon spin density matrix is taken 

in the form ege*^ — kl^k!^/\kT\'^, where kr is the compo- 
nent of the gluon momentum perpendicular to the beam 
axis. In the coUinear limit, when fc^ — ^ 0, this expres- 
sion converges to the ordinary e^e*"^ = —g^'^/2, while in 
the case of off-shell gluons it contains an admixture of 
longitudinal polarization. 

All algebraic manipulations with Feynman diagrams 
have been done using the computer system form [l^. 
We have checked that in the coUinear limit we reproduce 
the results of Refs. [13, [Il|. The full Fortran code 
for the matrix element is available from the authors on 
request. This process is also available in the hadron level 
Monte Carlo generator CASCADE [23|. Numerical results 
shown in the next section have been obtained using the 
AQ gluon distribution from [23| . 

In addition to the above, we also consider the pseudo- 
diffractive gluon-gluon scattering subprocesses repre- 
sented by the diagrams of Fig. 2. Despite the latter 
are of formally higher order in a^, they contribute to 
the events with large rapidity difference between the two 
J/f/) mesons and in that region can take over the leading- 
order 'box' subprocess. Our processes differ from the true 
diffraction in the sense that there occurs color exchange, 
and so, the rapidity interval between the two J /ij^ 's may 
be filled up with lighter hadrons (thus showing no gap in 
the overall hadron density) . Among the variety of higher- 
order contributions, the pseudo-diffractive subprocesses 
mentioned here are of our special interest as they poten- 



tially can mimic the DPS mechanism having very similar 
kinematics. 

The evaluation of the one-gluon exchange diagrams 
.g(fci) + g{k2) J/ipipi) + J/V'fe) + giks) + 9{ki) is 
straightforward, but the number of diagrams is rather 
large. There are six possible gluon permutations in the 
upper quark loop and six permutations in the lower 
loop. Besides that, we have to consider interchanges be- 
tween the two initial or two final gluons: g{ki) O 
5(^3) ^ 9(^4), thus ending up with 144 possible combi- 
nations. Note that the matrix element is free from in- 
frared singularities. This is due to the specific property 
of the quark loop amplitude which vanishes when any 
of the three attached gluons becomes soft. These calcu- 
lations have also been performed in the fct-factorization 
approach as described above. 

The two gluon exchange mechanism has been previ- 
ously considered in Ref. flil, where it was reduced to 
the production of J/ip pairs in photon-photon collisions 
[2^ by recalculating the appropriate color factor. We ba- 
sically follow the same way in our present analysis, but 
use an updated gluon density 1265 . 

Let us concentrate on the elementary g+g^J /tp+J /tp 
subprocess first. Only 16 of the different possible Feyn- 
man diagrams survive in the high-energy limit, which 
seems to be a suitable approximation for the conditions 
in discussion. The corresponding amplitude can be cast 
into the impact- factor representation (25| : 



,,2 J{9\i^y\i;i^,(i)J{9\2- 

is I a K ■ 



>Vx, 



(1) 



and the cross section reads 



daigg^VV;s) 
dt 
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Here the subscripts A; denote the helicities of the glu- 
ons g and vector mesons V , and q is the transverse mo- 
mentum transfer, t —q^- The overall color struc- 
ture of the reaction is described by the factor Nc = 
{N'l-if /[im'^[N'^-\% where = 3. We kept ex- 
plicit an effective gluon mass no which is responsible for 
soft QCD effects 27] and plays the role of regularization 
parameter. However, in the present study we can safely 
set it to zero, and the amplitude remains finite as the 
impact factors J vanish when k ±q/2. See also the 
discussion in (28j . 

At small t, within the diffraction cone, the cross sec- 
tion is dominated by the s-channel helicity conserving 
amplitude. In this case, the explicit form of the impact 
factor is 

Jigx^V,;n,q)^Sx,rV^s ! "^^"'^^^t^'f dzd^k, 

J z[l - z){2ny 

(3) 

where ipi^z^ k) is the light-cone wave function of the vector 
meson and z is the light-cone momentum fraction carried 
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by the heavy quark. Neglecting the intrinsic motion of 
the quarks we set ^p{z, k) — C 6{z — i) 5'^'^\k), where the 
normalizing constant C is adjusted to the J/il} leptonic 
width and is related to the radial wave function at the 
origin as = 12t:^ / {N'^m^^)\n{Q)\'^ . Within the above 
approximation, we have 

1 4 



m3,/4 + m? 



(4) 



Including the quark intrinsic motion would decrease the 
amplitude; for a more detailed analysis see [1^. As it 
will become clear from the numerical results, the non- 
relativistic approximation is sufficiently accurate for our 
purposes. 

The cross section for the two-gluon exchange contri- 
bution to the p + p — >■ J/V' + J/V" + X reaction (see 
Fig. 2) is calculated in the collinear approximation with 
MSTW2008(NLO) gluon distribution function [2^ and 
the factorization scale /i^ = ml, where rrit is the J/^ 
transverse mass. The elementary g+g — > J /il>+J /ip cross 
section can be easily calculated in the high-energy ap- 
proximation similarly to how it was done for the 7 -I- 7 ^■ 
J / -\- J / reaction [28|. The corresponding cross section 
is proportional to as(/x^), and therefore depends strongly 
on the choice of the renormalization scale. In the calcula- 
tion presented here we take = . In the high-energy 
approximation, the matrix element is merely a function 
of the transverse momentum q of one of the J /'ip 's. This 
cannot be true at low subprocess energy, close to the J /tp 
J /ip threshold. Here one must take into account also the 
longitudinal momentum transfer. We therefore replace 
q^ by the exact t or u for the t and u diagrams respec- 
tively. We neglect here the possible interference between 
the box diagram and the two-gluon exchange mechanism, 
which is formally of lower order than the square of the 
two-gluon amplitude. However, firstly, the addition of 
box and two-gluon exchange amplitudes is not warranted 
without the consistent evaluation of ag-corrections to the 
box. Secondly, it will become obvious from the numeri- 
cal results, that the two-gluon mechanism is exceedingly 
small in the region of invariant masses dominated by the 
box mechanism. 



B. DPS contributions 

Under the hypothesis of having two independent hard 
partonic subprocesses A and _B in a single pp collision, 
and under further assumption that the longitudinal and 
transverse components of generalized parton distribu- 
tions factorize from each other, the inclusive DPS cross 
section reads (for details see, e.g., the recent review [2| 
with many references to prior works listed therein) 



C^DPS 



Ceff 



A B 

m o-spsO"sPS 

2 CTcff 



dH{T{h)y 



(5) 



(6) 



where T(b) = J /(bi)/(bi— b) d^&i is the overlap func- 
tion that characterizes the transverse area occupied by 
the interacting partons, and /(b) is supposed to be a uni- 
versal function of the impact parameter b for all kinds of 
partons with its normalization fixed as 



/(bi)/(bi-b) dHi dH - / T(b) d^fc ^ 1. (7) 



The inclusive SPS cross sections tig^pg and a^pg for the in- 
dividual partonic subrocesses A and B can be calculated 
in a usual way using the single parton distribution func- 
tions. The symmetry factor m equals to 1 for identical 
subprocesses and 2 for the differing ones. 

These simplifying factorization assumptions, though 
rather customary in the literature and quite convenient 
from the computational point of view, are not sufficiently 
justified and are currently under revision jlj. Neverthe- 
less, we restrict ourselves to this simple form ([S]) regard- 
ing it as the first estimate for the DPS contribution. The 
presence of correlation term in the two-parton distribu- 
tions results in reduction [29l - l3]| of the effective cross 
section (Teff with the growth of the hard scale, while the 
dependence of (Teff on the total energy at a fixed scale 
is rather weak fsO]. Thus, in fact, we obtain the lower 
bound estimate for the contribution under consideration. 
The CDF 32„ ^] and DO J4j] measurements give tToff — 
15 mb, that constitutes roughly 20% of the total (elastic 
+ inelastic) pp cross section at the Tevatron energy. We 
will use this value in our further analysis. 



III. RESULTS AND DISCUSSION 

We start with discussing the role of kinematic restric- 
tions on the J/'0 transverse momentum. Shown in Fig. 
3 are the fractions of SPS events surviving after impos- 
ing cuts on priip)- Dashed line corresponds to requiring 
PTi'4')>PT,min for Only One (arbitrarily chosen) J/ip me- 
son with no restrictions on the other. Were the two J/V^'s 
produced independently (like in the DPS mode), the 
probability of having PT{'>P)>PT,min for the both J/ip's 
simultaneously could be obtained by just squaring the 
single-cut probability (dash-dotted curve in Fig. 3). On 
the contrary, in the naive the back-to-back kinematics, a 
cut applied to any of the two J/ip^s would automatically 
mean the same restriction on the other, thus making no 
effect on the overall probability (dashed curve). 

The explicit calculation (solid curve) lies in between. 
In the region pT,min < 4 GeV the solid and dash-dotted 
curves almost coincide, thus showing that the two J/ip's 
are nearly idependent. With stronger cuts on px^tp), the 
curves diverge showing that the production of J/ip's be- 
comes correlated. 

Another illustration of this property is given by the dis- 
tributions in the azimuthal angle difference da{iptp)/dAip 
exhibited in Fig. 4. The distribution looks flat for the 
unrestricted phase space (the upper plot), but tends to 
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concentrate around /S.ip ~ tt when the cuts on pt{iP) be- 
come tighter (the middle and the lower plots in Fig. 4.) 
In principle, one could get rid of the SPS contribution by 
imposing cuts like pritp) > 6 GeV, A(fi < tt/A, but the 
DPS cross section would then fall from tens of nanobarns 
to few picobarns. We can conclude that the SPS and DPS 
modes are potentially distinguishable at sufficiently high 
Pt{iP): but the production rates fall dramatically, and 
so, the practical discrimination of the production mech- 
anisms remains problematic. 

Now we turn to rapidity correlations explained in Figs. 
5 and 6. In the case of independent production (the 
DPS mode), the distribution over Ay is rather flat (dash- 
dotted curve in Fig. 5), while in the case of SPS 'box' 
contribution (dotted curve in Fig. 5) it is concentrated 
around Ay ~ and does not extend beyond the interval 
\Ay\ < 2. 

In Fig. 5 we also show pseudo-diffractive contribu- 
tions from the one- and two-gluon exchange processes of 
Fig. 2. As it was expected, these processes lead to rela- 
tively large Ay and even show maxima at Ay ~ ±2. This 
corresponds to a typical situation with one J/ip moving 
forward and the other one moving backward, following 
the directions of the initial gluons. The minimum for 
the two-gluon exchange g + g ^ J/ip + J/i' subproccss 
is a consequence of the educated guess correction of the 
high-energy formula at low energies as discussed above. 

At the same time, the absolute size of the one-gluon 
exchange cross section is found to be remarkably small. 
There are several reasons taking credit for this small- 
ness. First, is the presence of two extra powers of as- 
Second, is just the large typical rapidity difference that 
makes the invariant mass of the final state relatively 
large: M^^{Ay=2) /M^^{Ay=0) ~ cosh(Ay/2). In turn, 
larger masses mean larger values of the probed x, and, 
accordingly, lower values of the gluon densities. 

The third and the most important reason lies in the 
color factors. The color amplitude of the first dia- 
gram in the first row of Fig. 1 reads tr{T"T''T''T^} = 
[{Nl - l)/(4iVc)](5'^'' (where T stand for the SU{2,) gen- 
erators). After taking square and summing over the 
initial gluon colors a and b it gives [It^"**]^ — 32/9. 
Similarly, for the diagrams in the second row we have 
[tr{T''T''T'^}f'''^f = = = 9/2 and 

For comparison, the color amplitude of the first dia- 
gram in Fig. 2 reads jd"-'''^' ■ ^d!"''^, and the terms con- 
taining 1/°'^'^ or ■jf'"^'^ disappear because of cancellations 
between the different diagrams. This yields after squar- 
ing 



(jV,^-l)(A^,^-4)^ 

256 m 



1 200 
256 ~9~ 



0.1. 



tive): 

[d"" d'' 



'd'"" d 



h ]bce-\ 



(7\^-1)(A7-4)(A;--12) 



The color interference term is even smaller (and ncga 



^ _ -20 

256 ^ 512 iV2 ^ 256-3' 

A remark for those loving the 1/A^c expansion: all the 
considered contributions are of the same order in Nc- 

The same suppression factors apply to the two-gluon 
exchange as well, but there is yet another suppressing 
mechanism specific for the one-gluon exchange process. 
It comes from the fact that the amplitude vanishes when 
any of the final state gluons becomes soft (this property 
makes the process infrared-safe, as we have mentioned al- 
ready). Recall that by the same token the inclusive pro- 
duction rates of J/tp and Xc mesons become comparable 
to each other in spite of the hierarchy of the wave func- 
tions (5- wave compared to P-wave). As a consequence, 
although the two-gluon exchange g+g J/t/j+J/ip and 
one-gluon exchange g+g J/tp+J/ip+g+g processes 
are of the same QCD order, their magnitudes are con- 
siderably different. 



IV. CONCLUSIONS 

We have considered the production of J/ tp pairs at the 
LHC energies via SPS and DPS processes taking into 
account several possible contributing subprocesses. We 
find it rather difficult to disentangle the SPS and DPS 
modes on the basis of aziniuthal or transverse momentum 
correlations: the difference becomes only visible at suffi- 
ciently high pt, where the production rates are, indeed, 
very small. 

Selecting large rapidity difference events looks more 
promising. The leading order SPS contribution is local- 
ized inside the interval |Ay| < 2 (and continues to fall 
down steeply with increasing |Ay|), while the higher or- 
der contributions extending beyond these limits are heav- 
ily suppressed by the color algebra and do not constitute 
significant background for the DPS production. 
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FIG. 1: Examples of Feynman diagrams representing the 
leading-order gluon-gluon fusion subprocess gg J /^J /^. 
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FIG. 2: Examples of Feynman diagrams representing the pro- 
duction of J/-i/) pairs in pseudo-diffractive gluon-gluon scat- 
tering via one-gluon exchange and two-gluon exchange mech- 
anisms. 
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FIG. 3: Fraction of the production cross section left after 
imposing cuts on the J/ij} transverse momentum, see explar- 
nations in the text. 




FIG. 4: Azimuthal angle difference distributions after impos- 
ing cuts on the J/xp transverse momenta. Upper plot, no cuts; 
middle plot, pT,min = 6 GeV; lower plot, pT,min = 12 GeV. 
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FIG. 5: Distribution over the rapidity difference between 
,J/ip mesons. Dotted curve, SPS 'box' contribution; dashed 
curve, one-gluon exchange contribution multiplied by 1000; 
solid curve, two-gluon exchange contribution multiplied by 
25; dash-dotted curve, DPS production. 




FIG. 6: Double differential distribution da/dy{ijji)dy{i}!2) for 
the SPS production mode. 



